HIGHER CORRELATIONS OF DIVISOR SUMS RELATED TO 
PRIMES III: /(-CORRELATIONS 
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Abstract. We obtain the general fc-correlations for a short divisor sum re- 
lated to primes. 



1. Introduction 

This is the third in a series of papers concerned with the calculation and ap- 
plications of higher correlations of short divisor sums that approximate the von 
Mangoldt function A(n). We return to the short divisor sum approximation of the 
first paper defined for R > 1 by 

(1.1) A H (n) = £>(d)kg(iZ/d). 

d\n 
d<R 

The correlations we are interested in evaluating are 

N 

(1.2) S k (N, j,a) = J2 A n( n + nT 1 A u{n + hT 2 • ■ ■ A R (n + >)"' 

n=l 

and 

N 

(1.3) S k {N,j,a) = A fl (n + A R {n + 32 ) a - ■ ■ ■ A R (n + j r _ 1 )"-'A(n + j r ) 

n=l 

where j = (ji, J2, ■ ■ ■ , jr) a nd a = (ax, a>2, ■ ■ ■ a r ), the ji's are distinct integers, 
aj > 1 and ^2i—i ca = k. In (1.3) we assume that r > 2 and take a r — 1. We also 
define 

iV 

(1.4) S 1 (N,j,a) = Y. A ( n + ^~ N 

n=l 

if |ii I < AT by the prime number theorem. In the first paper in this series we 
evaluated Sk(N,j,a) and Sk(N,j,a) for 1 < k < 3, and also showed how for any 
k the two correlations are related. In this paper we will consider the general case. 
For j = ■ ■ ■ ,j r )i where the ji's are distinct integers, we define the singular 

series 



(1.5) 



^nHH 1 -^) 
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where v P (j) is the number of distinct residue classes modulo p that the jVs occupy. 

Theorem 1.1. Given k > 1, let j = (ji, j2, ■ • ■ , jr) and a = {0-1,0.2, ■ ■ -a r ), where 
the ji 's are distinct integers, and a,i > 1 with Ylj—i o>i = k. Assume max^ \ji\ < R 
and R > 2. Then we have 

(1.6) S k (N,j,a) = (C k (a)6(j) + o k (l))N(\ogR) k - r + 0(R k ), 

where the C k (a) are constants that are computable rational numbers. For k < 4 we 
have 

Ci(l) = 1, 

C 2 (2) = 1, C 2 (l,l) = l, 

C 3 (3) = \, C 3 (2,l) = l, C 3 (l,l,l) = l, 

C 4 (4) - C 4 (3,l) = ^ C 4 (2,2) = l, C 4 (2,l,l) = l, C 4 (l, 1,1, 1) = 1. 

Denoting C k {k) as C k , all of the constants Cfc(a) with a = (ai, 0,2, . . . , a r ) are de- 
termined from the constants C k by the formula 

r 

(1-7) C fc (o) = JJC« 4 . 

i=l 

We actually prove a more precise result than (^^), but for applications the result 
stated here is sufficient. 

If we wish to apply Theorem [O] and obtain explicit numerical results we need 
to know the values of the constants C k - These constants arc defined in terms of 
multiple integrals which may be evaluated by residue calculations. These residue 
calculations becomes increasingly difficult as k gets larger. We will obtain in this 
paper the values 

C\ = 1, C2 = 1, C3 — -, C4 = -. 

David Farmer has computed further values. He obtained the value C4 given above 

and also has found 

11065 / , rTror „ 11460578803 „„ Wrt „ 
C 5 = "214- = -67535 . . . , C e = ^ = .66709 .... 

From the first paper in this series, we can evaluate the mixed correlations (1.3) 
as an immediate consequence of Theorem 1.1. We need the following form of the 
Bombieri- Vinogradov theorem. Let 

x 



(1.8) E(x;q,a)=tp(x;q,a)-[(a,q) = l]- 

where [-P(x)] is the Iverson notation where square brackets around a true-false 
statement takes the value 1 if the statement is true and if the statement is false. 
Suppose for some fixed < •& < 1 that 

(1.9) V max \E(x;q,a)\ < 

— a loff x 

l<q<x<>-' (a,q) = l U& 

holds for any e > 0, any A — A(e) > 0, and x sufficiently large. This is a weak- 
ened form of the Bombieri- Vinogradov theorem if d = |, and therefore (1.9) holds 
unconditionally if $ < h. Elliott and Halberstam conjectured (1.9) is true with 
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§ = 1. (The constant $ is often referred to as the level of distribution of primes in 
arithmetic progression.) 

Theorem 1.2. Given k>2, let j = (ji,j2, ■ • ■ ,jr) an d a = (&i)&2 5 ■ • -a r ), where 
r > 2, a r = 1, and where the ji 's are distinct integers, ai > 1, and 53j=i a » = 
Assume max, < i? e and i? > A^ e . T/ien we have, /or iV e < i? < iV mm( ^'l?^ 
where (|l.9|) Zio/ds lozi/i 



(1.10) S k (N,j,a) = (C fc (a)e(j) +o fc (l))iV(logfl) fc - r . 

where Ck (a) are the same constants as in Theorem \1 . l\. 



We can also improve the range of i? in Theorem L2 to i? <C N 1 ^^ 11 by applying 
the methods of this paper directly to the mixed correlations, but this improvement 
is only significant in the conditional case when 1 — i < $ < 1. We will delay 
our applications to primes to the next paper in this series. To optimize some 
applications we will need to generalize slightly the previous results by allowing 
for different truncation lengths in the correlations. (We chose not to present the 
result initially in this generalized form to avoid confusing the reader with further 
notational difficulties which are easy to insert at a later point.) 



Our proof of Theorem 1.1 is largely self-contained and can replace the more 
complicated (but maybe conceptually easier) proofs given in the first paper in the 
case k = 2 and k — 3. We will first treat the case of pair correlation k = 2 in detail 
in the next section. The following three sections will handle the case k — 3 of triple 
correlation, where we give a more detailed proof then strictly necessary in order to 
provide examples for the general case. In the sixth section we prove the general 
case, and the following section discusses the computation of the constants Ck- In 
the final section we obtain the generalized result mentioned above for use in some 
applications. 

Notation. We retain the notation from our first paper which we will define 
here when it first appears. However, there are a few conventions we will use so 
frequently throughout the paper that they need to be m ent ioned immediately. We 
will use the Iverson notation which was defined below ([0]). We will take e to be 
any sufficiently small positive number whose value can be changed from equation 
to equation, and similarly C, c, and d will denote small fixed constants whose 
value may change from equation to equation. For a vector j = (ji,j2, ■ ■ • ,jk) we 
let Hj'II = max|ji|. A dash on a summation sign J^' indicates all the summation 
variables are relatively prime with each other, and further any sum without a lower 
bound on the summation variables will have the variables start with the value 1. 
Empty sums will have the value zero, and empty products will have the value 1. 
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2. Pair Correlation 



In this section we will prove the pair correlation case k = 2 in Theorem 1.1. The 
proof is more difficult than our earlier proof in (!]] , but the method is self-contained 
and generalizes nicely to higher correlations. The result we prove in this section is 
stated in the following theorem. 

Theorem 2.1. Let 

N 

(2.1) S 2 (j) = J2 A n(n)A R (n + j). 

n=l 

Then with T> a constant and c' a small positive constant, we have 

(2.2) 5 2 (0) ^ N\ogR + VN + 0{Ne~ c ' VT °^ R ) + 0(R 2 ), 
and, if j ^ 0, and letting j — (0,j), then for log \2j\ <C logR, 

(2.3) S 2 (j) = N&(j) + 0{Ne- c '^ n ) + 0(R 2 ). 



We have 



A' 



5 2 (i)=^A fl (n)A R (n + j) 

71 = 1 

= E E KdMe)log(R/d)\og(R/e) 



(2.4) n<N d,e<R 

d\n, e|n+j 

= Y, Kd)fi(e)\og(R/d)log(R/e) £ 1. 

d,e<R n<N 

d\n 
e\n+j 

The divisibility conditions d\n and e\n + j can only be satisfied if (d,e)\j, in which 
case n will run through a unique residue class modulo [d, e] . Hence 

(2.5) J2 l = [W e )ii](#r + °( i : 

n<N V ^' eJ 

d, | n 
e\n+j 

and therefore 

(2.6) S 2 (j) = N ]T ^f±\og(R/d)\og(R/e) + 0(R 2 ). 

d.e<R [a ' 6J 
{d,e)\j 

To handle the least common multiple in the sum above, we let d = b\b\ 2 and 
e = b 2 b\ 2 , where 612 = (d,e). Then [d, e] — b\b 2 bi 2 and b%, b 2 , and 612 are pairwise 
relatively prime. (This notation conforms with the notation we will use in the 
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general case.) Hence we have 

b\b 2 b\ 2 " b 61612 6 2 6i 



(2.7) b 2 b 12 <R 

= NT 2 (j)+0(R 2 ), 

where the dash on the summation sign indicates the summation variables are pair- 
wise relatively prime with each other. At this point in the earlier proof in 0] we 
summed over each variable individually. Here instead we apply the formula, for 
m > 2 and c > 0, 

■:<" 1:! /" ^ , f0, if0<*<l, 

1 J 27rf J c _ ioo 8 m \ (l0gx) m -\ ifx>l. 

Letting m = 2 in this formula, and denoting the vertical line contour c + it, —00 < 
f < 00 by (c), we see that for c\,c 2 > 0, 

(2.9) T 2 (j) = — / / F( Sl , S2 ) — —d Sl ds 2 , 

{2m) 2 J J si 2 s 2 2 

(C2) (d) 

where 

ju(6i)^(6 2 )At 2 (6i 2 ) 



(2.10) F( Sl ,s 2 )= J2' 



b l 1+Sl b 2 1+S2 b 1 ~ 1+Sl+S2 ' 



Letting Si = + zii and s 2 — a 2 + zt 2 , we see F(si, s 2 ) is analytic in si and s 2 for 
o"i > and 02 > 0. To analytically continue F(si, s 2 ) to the left we start with the 
product representation 

(2.11) ^ i _ ( 

^(S1,S2) = JJ f 1 - 1+Sl - p 1+S2 j ]^[ f 1 - pl+s~i ~ pl + s 2 + pl+si+ S 2 J ■ 

vh ' p\j 

There are two cases to consider. If j — then the first product is empty and 
F(si, s 2 ) is equal to the second product over all primes, while if j ^ then the sec- 
ond product is a finite product, and the analytic continuation of F(si, s 2 ) depends 
on the first product. We will deal with each case separately. 
Case 1. Assume j = 0. In this case, we have 

(2.12) F( Sl ,s 2 ) = l[(l ' ' ' 

p ^ 

Since for Re(s) > 1 

(2-i3) c( S )=n( i -i 

and, for min(l + 01, 1 + <j 2 , 1 + o\ + a 2 ) > 0, 
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we see that 

(2-14) F(s u s 2 ) = 77TX ± WxM Sl ' S2) 

where 

(2.15) h(s 1 ,s 2 )='[[ 

p 

and 



1 



v 



1 1 i 1 



pl+Sl+S2 / ' 



Msi, S2 )=n(i+o(- 2 V)) 

with v = (o"i + era ) + min(0, eri, (72, <7i + 02). Hence we see h(s%, s 2 ) is analytic in s\ 



and S2 for <ri > — \ and cr 2 > —j. Therefore ( 2.14J ) provides the analytic continu- 
ation of F(s\,s 2 ) into this region and shows that |/i(si , S2) | < 1 as |ti|,|t2| ~~ * °o- 
We also see that h(0,0) = 1. 

It should be kept in mind that while F(si,s 2 ) is a function of two complex 



variables, we will at each stage in our calculation of (2.9) always fix one variable 
and treat F as a function of the other single complex variable. Thus, in the analysis 
that follows we assume in each equation that one of the complex variables s± or s 2 
is fixed. When we treat the general ^-correlations we will always fix all except one 
complex variable in F at each stage in the calculation. 

We now isolate the dominant part of F(s\, S2) around si = 0, s 2 = 0. Since 

00 

( S - i)c(s) = 1 + y> nS " 



71 = 1 



, {si + s 2 )C{l + Sx + s 2 ) , A S 1 S 2 
F(s 1 ,s 2 ) = — — — — — — -h(si,s 2 ) 1 



is analytic in the entire s-plane, we write 
(2.16) 



siC(l + si)s 2 C(l + Si) ' Js!+s 2 
sis 2 



= k(si,s 2 ) 

si + s 2 

We see 

(2.17) fc(si,0) = fc(0,s 2 )=l, 

and fc(si,S2) is analytic everywhere in si and s 2 except for poles at the zeros of 
C(l + si) and C(l + s 2)- To avoid these poles, we use a classical zero-free region 
result. By Theorem 3.11 and (3.14.8) of || there exists a small positive constant 
C such that £(<r + it) 7^ in the region 

(2-18) cr>l-- — ,. C . r 
V ' ~ Iog(|*|+2) 

for all t, and further 

(2.19) ((<j + it) L— « Iog(|*| + 2), - 1 « log(|t| + 2) 

(7 — 1 + it £(<7 + it) 

in this region. (There are stronger results but this suffices for our needs.) Hence 
k(si, s 2 ) is analytic in si and s 2 in the region 

C C 

(2.20) 0-1 > - ct 2 > 



Iogdtil+2)' z " Iog(|t a |+2)' 
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and we also see by ( 2. IE ) that in this region 
(2.21) F(si,s 2 ) < log(|*i| + 2)log(|i 2 | +2) max i 



si + s 2 



,log(|ti+t 2 | + 2)). 



Returning to ( p.9|) in this case, we have 
1 



(2.22) 



r a (o) 



(2ttz) 2 



k(si, s 2 )- 



■ <isi ds 2 , 



sis 2 (si + s 2 ) 

(02) (Cl) 

where we now take c 2 > _£i> 0. To evaluate this integral, we move the contours to 
the edge of the region ( 2.20 ) and evaluate the residues. Let C denote the contour 
given by 

C 



(2.23) s = - - — — — — + it. 

log(|*|+2) 

and let Cj denote the same contour except with C — Cj. Since we want to avoid 
integrating over the pole si = —s 2 of the integrand on the contour we move to, we 
will move (a,-) to Cj and choose C 2 = \C\ so that the contours C\ and C 2 are well 
separated with 

, 1 

log(|si| + 2) 

on these contours. Now, provided \s\\ 3> 1, |s 2 | S> 1, and \s\ + s 2 | 3> i og (| s L 1 | +2 ) , 



(2.24) 



F( Sl ,s 2 ) 



R 



S1+S2 



Si 2 S 2 2 



< R° 



r2 log 2 (|t 1 | + 2)log 2 (|i 2 |+2) 



(1^1 + 2)2(^1+2)2 • 

in the region ( [2.20 ). This bound shows the integral in ( 2.22 ) is absolutely conver- 
gent, and, since we can easily arrange (ci) and (c 2 ) to meet the conditions necessary 
for this bound, which are also satisfied on C\ and C 2 , justifies the residue calcu- 
lations that follow. Moving first (ci) to C%, we encounter a simple pole at s% = 
with residue by (2.17), and hence we obtain using (2. 



T 2 (0) = 



1 

27ri 

lop 



R s ? 



1 



(2m) 



i\2 



k(si,s 2 ) 



sis 2 {si + s 2 ) 



dsi ds 2 



(c 2 ) (c 2 ) Cx 

1 f f R Sl+S2 

to -\2 / / k (si,s 2 ) -. ■ rds 1 ds 2 . 

(2ttiY J J sis 2 (si + s 2 ) 

(C2) Cl 

In the double integral above we move the contour (c 2 ) to the left to C 2 and encounter 
two simple poles at s 2 — and s 2 = — Si, with residues 

R S1 , fc( Sl ,- Sl ) 



and 



respectively. Hence 



1 



Si" 



(2.25) T 2 (0) = logi?+2?+ — / — dsi 



2ni 



R s 



where 
(2.26) 



1 

2tt7 



(27Ti) 



F(si, s 2 ) — — Tj, dsi ds 2l 



Si 2 S 2 2 



c 2 d 



k(si, -si) 



ds\. 
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The first integral in (2.25) is zero since we may move the contour to the left without 
encountering any singularities to — oo where it vanishes. Using (2.24) the remaining 
integral in (2.25|) is bounded by 



c log 2 (|t| + 2) , 

r io 8( iti+2) & y l Z dt 



t 2 + i 



The integral here is, for any w > 2, 



< log w I R l °sd'i+ 3 ) dt + 



t 2 



dt -C u>(log 2 w)e~ 



= logH fog W 



and on choosing logu; = ^clogR this is 

< (clogi?)e^ Vclog " 



< e 



We note for future use that the above argument shows that with s = a + it, 
and B a constant 



w - io g (|t|+2; 

(2.27) 
since here 

We conclude that 
(2.28) 



^ iogB ; |s i +2) ^«e-™ 



J_ <<; log(| S | + 2) 



2+ t : 



T 2 (0) = log R + V + 0(e 



which by (2.7) proves the first part of Theorem 2.1 . 

For computation purposes we note that the main term in the above analysis is 
obtained directly from 



(2.29) 



1 



(2ttz) 2 



R S1+S 



■ dsi ds2, 



SlS 2 (Sl + s 2 ) 

(ca) (ci) 

since the main term in ( [2.22 ) occurs from the residues with fc(0,0) as a factor. Wc 
make the change of variables Wi = (logi?)si and w% — (log R)s2- Since we can 
take Cj/ log R in place of Cj for j = 1, 2 in the double integral, we conclude that the 
main term is 

(2.30) C 2 (2)logi?, 
where, on returning to the variables si and s%, 



(2.31) 



C 2 (2) 



1 



(2ni) 



= Sl+S2 



(C2) (Cl) 



SlS2(si + ,S 2 ) 



(isi C?S2- 



We will discuss this integral and the computation of Cfc(fe) in Section 7. 

Case 12. j =/= 0. In this case, the first product in (2.11) is the dominant one, and 
therefore 

(2.32) F(s u s 2 )= ai + s J ai + s2) h jM 
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where 
(2.33) 

"'^ - S (d-^)(i-^)J n ( (i-^xi-^) j ■ 

and we see hj(s\, s 2 ) is analytic in the region min(l + a\ + 02, 1 + 01, 1 + 02) > 0, 
or in particular if ci, (T2 > — ^. Hence F(s±, s 2 ) is analytic in the region ( ^20| ). We 
also see 

p/fj V 1 P j P \j v P 

where, since j ^ 0, v p (j) is 1 if p|j and = 2 if p /j. Thus f p (j) = v p {j) when 
j = (0, j), and we have shown that 

(2.34) ^(0,0) = 6(j). 

We now move the contour integrals in (|2.9| ) to the left as before and encounter 
simple poles at S\ = and s 2 = which give 

(2.35) T 2 (j) = + -— j / F(s 1 ,s 2 )^ r ^ds 1 ds 2 . 

c 2 £1 

By the estimate preceding (3.5) of Q we have for <5 > 

(2.36) [] f 1 + °b^f )) « cx P(^( l0 S \ 2 J\) S ) 

for some positive constant D, and therefore in the — \ < a\, 02 < we see that 

^(^i^2)«n( 1 + Q ( p i + . 1 1+CT2 )) «exp(£>(log|2j|)-— 2 ). 
p\j 

We conclude by ( ^.19| ) and ( |2.27 ) that the double integral in (2.35) is 

« (exptDOogM)^) £ j^il2 ^e-'^ 
on taking (as we may) Ci +C2 < | log 2 and using log \2j\ <§; logi?. This completes 



the proof of Theorem 2.1 



3. Initial treatment of triple correlation 

We now consider the triple correlations of A#(n). The procedure is the same as 
for the pair correlations only the calculations are more complicated. We have been 
more detailed then necessary in order to provide motivation for the general case. 
We will prove the following theorem. 

Theorem 3.1. For constants T>q and T>\ we have 
(3.1) 

N 3 

^(O) = ^A^(n) 3 = -AHog 2 R + V X N log R + V Q N + 0(Ne~ c VT °^ Jl ) + 0(R 3 ). 

n=l 

For j = {31,32), h ¥=h 7^ 0, with log2||j|| < logi?, and a = (2,1) we have 
(3.2) S 3 (N,j,a) = 6(j)N\ogR + V 2 (j)N + 0(Ne- c '^^) +0(R 3 ), 
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where T> 2 {j) is an arithmetic function which satisfies the boundT>2(j) <C exp(||j|| 5 ) 
for any fixed S > 0. For j = (ji, j 2 , 33), ji ^ h ^ h ^ with log2||j|| < logR, 
and a = (1, 1, 1), we have 



(3.3) 
Let 

(3.4) 



S 3 (N,j, a) = &(j)N + 0(7V e -c'v^) + 0( ^ R 3y 



N 



S3 (ii , 32 , 33 ) = Ar ( n + - ?1 ) Ar ( n + - ?2 ) A « ( n + is) • 



n=l 



Expanding, we have 

S3{h,h,h)= v(d 1 )\og{R/d 1 )ii(d 2 )\og{R/d 2 )ii(d3)\og{R/d3) J2 L 

di,d2,d3<R 



n<N 
n+ji 
n+j 2 
n+j 3 



The sum over n is zero unless (di, d 2 )\j 2 — ji, (di, ^3)^3 — ji, and (d 2 , ^3)^3 — j 2 , 
in which case the sum runs through a unique residue class modulo [di, d 2 ,d3], and 
we have 

£ 1 v 



n<N 
n+ji 
n+j 2 
n+33 



[di,d 2 ,d 3 ] 



0(1). 



We conclude, letting 

(3.5) 

that 
(3.6) 



Li(R) - log 



d, ' 



E 



di,d 2l d3<R 
(di,d 2 )|j2— Ji 



[di,d 2 ,d 3 ] 



L 1 (R)L 2 (R)L 3 (R) + 0(R 3 ) 



(di 
(d 2 



<i3)|j3-jl 
d.3)b'3-j2 



We now decompose di, d 2 , and d3 into relatively prime factors 

dl = &1&12&13&123, 

d 2 = b 2 bi 2 b 2 3bi23, 

d3 = 63 &13&23 ^123 j 

where 6 X is a divisor of the dj's where i occurs in the digits of \- Since the dj's are 
squarefree, these new variables are pairwise relatively prime. We will let V denote 
the set of 6 x 's which satisfy the conditions 



&12&123|j2 - jl, 6136123IJ3 - jl, &23&123|.?3 ^ 32- 
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Then we have 
(3.7) 

MMM^M(^ 2 (MM 2 (MM 2 (ft23Hftl23) T /mT rz? w (m 
u v. C~t ^ D O162O3O12O13O23O123 

b3bl3623fcl23<fl 

V 



Now by (2.8) we have 

1 f f f R Sl R S2 R Ss 

(3.8) T a (j 1 ,j 2 ,ja) = , .nq / / / ^1,32, s 3 ) —J — 2 — 2-rfsirfs 2 ds 3 , 

(2tti) j J J J si 2 s 2 l s 3 

(ca) (C2) (ci) 

where 
(3.9) 

* {Si,S2,S3) ^ , 1+sii l+s 2 . l+s 3 . I+S1+S2L I+S1+S3L l+«a+*3l. l+si + s 2 +s 3 ' 
jj °1 02 03 012 °13 O23 O123 

We define J12, J13, and J23 to be the largest squarefree divisor of j 2 — ji, 33 — ji, 
and j 3 — j 2 respectively. Let 

(3.10) A 3 = (j 2 -ii)(j3 - ji){j3 ~ h), K= (J 12 ,J 13 ,J 23 ). 

Since j 3 — J2 = (is — ji) — (j2 — ji), we see that if a prime p /re, then it can divide 
at most one of numbers J12, J13, J23- Here for convenience we will extend the usual 
definition of the gcd by (a, 0) = a, and (0, 0) = 0. With s n = <J n + it n , 1 < n < 3, 
we have for <r„ > with 1 < n < 3 the product representation 

(3.11) F(s 1 ,s 2 , s 3 ) = hi(n)h 2 {K, Ji 2 )h 2 (n, Ji 3 )h 2 (n, J 23 )h 3 (A 3 ), 

where 
(3.12) 

hi{n) = h 1 (K,s 1 ,s 2 ,s 3 ) = 

nO 1 1 



pi + Sl pl + S 2 J5I + S3 pl + S!+S 2 pl+Sl+8 3 pl + S2+S3 pl + Sl+S 2 +S 3 



(3.13) 

h 2 (n,J ij ) = h 2 (n,J ij ,s 1 ,s 2 ,s 3 ) = Y[ 1 



P| Ji; 



pl + Sl pl + «2 J5I + S3 pl + Si+Sj 



for 1 < i < j < 3, and 

(3.14) /i3(A 3 ) = /i3(A3 ) s 1 ,S2,S3)= II f 1 -^? 



P/A3 



Jjl + Sl pl + S2 pl + S3 



There are three cases to consider, when j\ = j 2 = j 3l when j± = j 2 ^ j 3 (or a 
permutation of this) , and when ji 7^ j 2 j 3 . These cases correspond to the three 
numbers j 2 — ji, j 3 — ji, and j 3 — j 2 falling into 1, 2, or 3 different residue classes 
for all sufficiently large primes. We consider each case separately. 
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4. 5 3 (0) 

We consider the first case when j\ = j% = j 3 , where ( 3.11 ) becomes 
(4.1) 

^(81,82, S3) = /ll(0, Si, S 2 , S3) = 

1111 1 

V ' 



T](i-- - 

11 \ nl+si n l+ 



j)J-TS2 pl + «3 pi + Sl+S 2 pl + Sl+S 3 pl + S2+S3 pl + Sl+S 2 +S 3 



and we see 
(4.2) F(si,s 2 ,s 3 ) 



C(l + si + s 2 )C(l + si + s 3 )C(l + s 2 + s 3 ) 



C(l + Si)C(l + s 2 )C(l + s 3 )C(l + 81 + s 2 + s 3 
where h(0, 0,0) = 1, and 



■h(si,S2,s 3 ), 



(4.3) 



M S i, S2 , 53 )=n( i +°(^7)) 



for g\ + (T2 + 0-3 > — I , where, assuming eri < (T 2 < 03 , 

1/ < di + cr 2 + min(0, 01, <j\ + <7 2 ,<Ti + a 2 + (T 3 ). 

Hence h(si, s 2 , s 3 ) is analytic in s\, s 2 , an d S3 for 01, cr 2 , cr 3 > — ~, and \h(si, s 2 , s 3 )| <C 
1 as 1^1, \t 2 \, \t 3 \ — > 00. Next, as in ( 2.16 ) we define fc(si, s 2 , s 3 ) by 

sis 2 s 3 (si + s 2 + s 3 ) 



(4.4) 

Hence 
(4.5) 



-F(si,s 2 ,s 3 ) = /c(si,s 2 ,s 3 ) 



(si + s 2 )(si + s 3 )(s 2 + s 3 ) 



fc(si,s 2 ,0) = fc(si,0,s 3 ) = fe(0,s 2 ,s 3 ) = 1 
and as before k(s\, s 2 , S3) is analytic in the region 

C 



(4.6) 



aj - log(|i j |+2)' 



1 < 3 < 3, 



log(N + 2)) 



for a small enough constant C, and in this region by ( |2.19 ) 
(4.7) 

f(si,s 2 ,s 3 )« (nio g 2 (w+2) ) n max d ]_ r . 

Vii / Vi<r<,< 3 i Si + ^i 

As before we let £7 denote the contour from ( [2.23 ) with the constant C = Cj. 
Starting with a constant C\ for which (2.18) holds, we take the constants Cj = 
hCj-i for the contour Cj, j — 2,3, which arranges the contours so that Cj is to 
the left of Cj+i and these contours are well-spaced from each other. We let 
(4.8) 

<T( \ 1 f [ f U ^ 1+S2+S3 ( Sl + S2 + S3) 

J(C 1 ,C 2 ,C 3 ) = / / / fc(Sl,S 2 ,S 3 ) ? ; ; - ; -ds 1 ds 2 ds 3 . 

{2-kiY J J J sis 2 s 3 (si 



(c 3 ) (c 2 ) (ci) 



i(si -I- s 2 )(si + s 3 )(s 2 + S3) 



and see by (3.8) that for c 3 > c 2 > c x > (or we can just take Cj = j) 
(4.9) T 3 (0) = J( Cl ,c 2 ,c 3 ). 
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We first move ci to the left to L\, and encounter only the simple pole at s\ = 
since the other two poles of the integrand at si = — S2 and s\ = —S3 are to the left 
of C\. By (|4.5D the residue of this pole at s\ — is 



and therefore by (|2.8|) 



s 2 2 s 3 2 



T 3 (0) = log 2 R + J(C U C2,c 3 ). 

Moving next (02) to the left to C2, we encounter simple poles as s 2 — and S2 = — Si 
(but not S2 — —S3) with residues 

R S1+S3 and k( Sl ,- Sl ,s a )R s * 



si 2 s 3 2 ' si 2 (si + s 3 )(-si + s 3 ) ' 

respectively, and hence 

1 f f R Sl+S3 

T 3 (0) =l g 2 i?+— — / / ^^ds 3 dsi 

(2tTi) 2 J J S! 2 S 3 2 

(4.10) £l (C3) , 

1 /• / fcs 1 ,-s 1 ,s 3 ^ 

(2my J J s 1 2 (s 1 + s 3 )(-si + s 3 ) 
£1 (c 3 ) 

The first integral above is zero since we can move C\ to the left to —00 where it 
vanishes without encountering any singularities. In the remaining terms we move 
(03) to £3; the second integral encounters a simple pole at S3 = — s\ but not the 
pole at S3 — si due to how we arranged the contours Ci, and in J we pass simple 
poles S3 = 0, S3 = — si, and S3 = — S2. We conclude 



(4.11) 

= lno- 2 R A- 

2m J 2si 3 



T 3 (0)=log 2 i? + -L / ^'-^^V 



£1 



k(si,- Sl , S3 )R s * If , 

ds 3 as 1 + /n _, N9 / / — 2 — 2" 2 1 



(2wi) 2 J J si 2 (s 1 +s 3 )(-s 1 + s 3 ) ° (2m) 2 J 7 si 2 s 2 

£1 £3 £1 £2 



s 2 , -s^i? 52 fc(si, s 2 , ~s 2 )R s 



(2m) 2 J J \si 2 (si + s 2 )(-si + s 2 ) s 2 2 (si + s 2 )(s x - s 2 ) 
£1 £2 

i7 (£1, £2, C 3 ) 



ds 2 dsi 



, 2 D 1 1 /" ^(si,-si,-si)i?^ Sl 

= l0g R +^~ / dSi - Ji +i 2 -J 3 +^(£ 1 ,>C2,£3)- 

£1 

The integral I 2 is zero since we may move the contours to the left to —00 where 
the integral vanishes. To estimate the remaining integrals, we note for Sj on Ci we 

have N > iog(k|+a) ' and 



|Sl±Sjl>> log(|s l | + 2)' W 
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Thus by @ and ( ^27|) 

J{L X ,L 2 ,L 3 ) « ^jf 



(4.12) 



« e 



Similarly we have in I\ 

k( Sl ,- Sl ,s 3 )R S3 = C(i + si + s 3 )C(i - gi + s 3 )R S3 

si 2 (si + s 3 )(-si + s 3 ) Sl 4 S3 2 C(l + si)C(l ~ si)C 2 (l + *s) 

^^log 4 (|t 1 | + 2)log 4 (|i 3 | + 2) 



h(si,-si,-s 3 ) 



i*in«3i 
r^ 



log 6 (|t 1 | + 2)log 4 (|i 3 | + 2) 



and thus (2.27) applies, and the same argument holds for I 3 . Finally, for the first 
integral in (4.11) we move the contour back to the right across the triple pole at 
si = to the contour L\* given by s\ — i og (| ( ^' 1 1 | + 2) + it which is the reflection of C\ 
over the imaginary axis. Thus this integral is 

'fc(si, -si, -si)Rr Sl 



(4.13) 



res 

si=0 



2 Sl 3 



1 

2ttI 



k(si,-s 1 ,-s 1 )R Sl 
2si 3 



The integral here can now be estimated as before by (2.27) (because 1 / 1 s 1 1 3 <C 
log(|ti| + 2)/|si| 2 on the contour) with the same error term. Finally, the triple pole 
has residue 



(4.14) res ( i?" 

si=0 



Si C(l - 2si)h(si,-si,-si 
Sl 6 C(l + si)C 3 (l-si) 



1 



•log 2 i? + X>i \ogR + T> , 



for constants T>q and T> x . Combining our results, we conclude 
(4.15) ' 1 



T 3 (0) = ~ log 2 R + V 1 logR + V + 0(e- c ^*). 



This completes the proof of (p.lf) in Theorem 3T 



For computational purp oses we note that the main term is obtained from terms 
with fc(0, 0, 0) which from (O) is 



1 



(c 3 ) (c 2 ) (ci) 



(Sl + S 2 + 5 3 )^ 1+S2+ ' 3 

sis 2 s 3 {s 1 + s 2 )(si + s 3 )(s 2 + s 3 ) 



dsi ds 2 ds 3 



On making the change of variables w\ = (logi?)si, w 2 = (log R)s 2 , w 3 — (log R)s 3 
and taking Cj / logi? in place of Cj for 1 < j < 3 this term becomes 



where, on replacing Wi by s 
1 



(4.16) C 3 (3) 



(27T4) 



(c 3 ) (c 2 ) (ci) 



C 3 (3)log 2 i?, 

(gi + s 2 + s 3 )e Sl + A2+S3 

SlS2S 3 (si + S 2 )(si + S 3 )(s 2 + S3) 



dsi ds 2 ds 3 . 



HIGHER CORRELATIONS OF DIVISOR SUMS III 



15 



5. S 3 (N 7 j,a) FOR a = (2, 1) and a = (1, 1, 1) 

We first consider the case ji = j 2 ^ j 3 . Let j =_ja — ji = ja — J2, and take J to 
be the largest squarefree divisor of j. Then in (3.1C) we have A3 = and K = J, 
and (3.11) becomes 

(5.1) 

1 1 



F( Sl ,S 2 , S 3 ) = h x (J, Sl ,s 2 , s 3 ) H ( 1 - — ^ - — j 
= hi(J, si,s 2 ,s 3 )h 2 (J, 0), 



i+s 3 



1+S1+S2 



where /ii is given by ( gl2| ) and h 2 

is given by ( |3 . 1 3| ) - On factoring out the dominant 
zeta factors from the product above we write 

C(l + 81 + 82) 



(5.2) F( Sl ,s 2 ,s 3 ) 



-hi(J, si,s 2 , s 3 )hj(si,s 2 , S3) 



C(l + *i)C(H-aa)C(l + «3) 
which by ( |5.l|) se rves to define hj(si, s 2 , s 3 ). By the same argument used for 
h(si,s 2 ) in (2.15) we see that hj(si, s 2 , s 3 ) converges absolutely and is analytic 
in the half-planes cr n > — c, 1 < n < 3, for c a positive constant (c = j is accept- 
able). We also see immediately that h\{J, 0, 0, 0)hj(0, 0, 0) = 6(j) , wit h j = (0, j) 
since we obtain the same product here as in the equation above ( 2.34 ). Next, we 
define 

(5.3) F(si,s 2 ,s 3 ) = k(si,S2,s 3 )- 



and let 
(5.4) 

^(ci,c 2 ,c 3 ) 



1 



(2m) 



(si + s 2 X(l + s 3 )' 
k(s llS2 ,s 3 )R^+ s - ( R s 



(c 3 ) (ca) (ci) 

so that by ( ]3~8| ) for c 3 > c 2 > Ci > 

(5.5) T 3 (j) = J(c u c 2 ,c 3 ). 



sis 2 (si+s 2 ) \s 3 2 ((l + s 3 ) 



ds\ ds 2 ds 3 



Now, in the region (4.6) for C appropriately small 

(5.6) hx(J, Si,s 2 ,s 3 )h j (s 1 ,S2, s 3 ) « [] ( 1 + °H=tf)) « cxp(^(log |2j|) 5 ), 



p\j 



by (2.36), and hence in this region provided in addition |sj| ^> 1, 1 < i < 3, and 

l S l + S 2| > iog(|A|+2) WeSee 

(5.7) F(s 1; s 2 , s 3 ) « e X p(D(log |2j|) 5 ) log 2 (|*i| + 2) log(|t 2 | + 2) Iog(|t 3 | + 2) 

We move the contour (03) to the left to C 3 which crosses the simple pole at 
S3 = and gives 



(5.8) 



1 



(2m) 



fc(si,s 2 ,0)i? Sl+S2 



ds 1 ds 2 +J(c 1 ,c 2 , C 3 ). 



sis 2 (si + s 2 ) 

(ca) (ci) 

The first integral may now be handled identically with how we handled the corre- 
sponding integral in (2.22) except the analytic factor k(si, s 2 , 0) has the value 

fc(ai,0,0) = fc(0,« a ,0) = 6(j). 
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We thus obtain 

(5.9) T 3 (j) = e(j)N log R + V (j)N + J (a, C2,Lz), 

where T>o(j) a lso sa tisfies the bound in (|]^). For the last term we proceed as we 
did following ( [2.24 ). We arrange in order the contours C\ and £ 2 to the right of 
C 3 . Moving first (ci) to C± we encounter a simple pole at si = 0, and next movimg 
(02) to the left to £2 we obtain residues at S2 = and S2 = — si- By (2.27), 



R S3 



S3 2 C(1 + S 3 ) 



dt 3 < e 



and therefore by (5.7) these terms satisfy the bound <C A^e c ^ lo s R for log \ 2j\ -C 
logR. The theorem now follows in this case from the estimate 



J(£i,£ 2 ,£ 3 ) « Ne 



— c'\/log -R 



which also follows from the bound in (5/7) and ( 2.27 ). 

We finally turn to the case when a = ( 1, f , 1 ). Thus we assume j\ ^ j 2 7^ J3, and 
in this case the only infinite product in (3.11) is h 3 (A 3 ). Removing the dominant 
zeta-factors we have 



(5.10) 
where 

(5.11) 



h 3 (A 3 ) = 



C(l + *l)C(l + *3)C(l + *3) 



h 4 (A 3 



^{5) = n 



><n(i 

p\S 



p 



1+si 



1 - 



P 



l+s 2 



1 - 



P 



I+S3 



C(l + Sl)C(l + S 2 )C(1 + S 3 ) 



hj(si,s 2 ,s a ), 



We conclude that 

(5.12) F( Sl ,s 2 ,s 3 ) 

where 

(5.13) hj(st, s 2l s 3 ) = hi(K)h 2 {K, Ji 2 )/i 2 (k, Ji 3 )/i 2 (k, J 23 )/i4(A 3 ) 



is analytic in the region (4.6). We evaluate T 3 in ( |3.8| ) by moving the three contours 
to the left edge of this region. The only singularities in this region are simple poles 
at si = 0, S2 — 0, and s 3 = 0, and therefore we obtain 



(5.14) 

T 3 (ji, j 2l js) = ^-(0,0,0) + 



JJS1+S2+S3 
F(sx, s 2 , s 3 ) — ^ — t, — - dsi ds 2 ds 3 



c c c 



Si 2 S 2 2 S 3 2 



^•(0,0,0) + 0(expp(log2||j||) 5 ))e 
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wh ere t he error term is obtained from the esti mate ( 2.19 ) together with the estimate 
in (5.6) which applies for hj{sx, S2, S3), and ( 2.27 ). Since 

w»)^H)nH)n(^)n('- 1,J 

p\ K V F/ pj(K V F/ p/A 3 M 1 p) 7 p\A 3 V 

(5.15) pI a 3 



n - 



^ 3 fi v j(p) 



where 



1 if J f\r 

Vj (p) = { 2, ifp|Aa,p ](k, 
3, ifp/Aa, 

we see that ^j(p) agrees with the definition of v p {j). We conclude by (3J3) that, 
for log 2H.7 1 1 «logi?, 

(5.16) S*(ju3*h) = &U) + 0(e- c '^^), 

which completes the proof of Theorem |3.lL 



6. The /c-correlations 
We start by considering the sum 



N 



(6.1) S(j) = J2^R(n + h )A R (n + j 2 ) ■ ■ ■ Afl(n + jk), 



n=l 



where j — (ji, j'2, . . . , jk) and the jVs are not necessarily distinct. Proceeding as 
before we find as in fl2.6|) and (3.6) that 



(6.2) .,*„«,<« [*.*■-,*] 

(<i r ,d 6 )|j s -jV,l<r<s<fe 

= NT k (j) + 0(R k ). 

We next decompose the ck's into relatively prime factors. Let V(k) be the set of 
all non-empty subsets of the set of k elements {1, 2, . . . , k} (This is just the power 
set with the empty set removed.) For B £ V{k) : we let Ts{k) denote the set of all 
members of V(k) for which B is a subset. Thus for example if k = 4 then 

7>{1, 2 }(4H{{1,2},{1,2,3},{1,2,4},{1,2,3,4}} 

Since the d^s are squarefree we now decompose them into the relatively prime 
factors 

(6.3) di= [] b m , l<i<k, 

meV {i} (k) 

where b m is the product of all the primes that precisely divide all the di 's for which 
i € m, and none of the other d's. This decomposition is unique and the 2 k — 1 
6 m 's are pairwise relatively prime with each other (Here we subscript with the set 
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m rather than a list of its elements, so that for example &12 in the earlier sections 
now is written ^{i^})- Letting T>(j) denote the divisibility conditions 



(6.4) 



n & » 



Jt - Js 



1 < s < t < k, 



we sec on substituting ( |6.3| ) into ( |6.2| ) and applying fl2.8| ) that, letting #m denote 
the number of elements of m, 



T k (j) 



e n 



(6.5) 

where 
(6.6) 

and 
(6.7) 



di,d 2 ,...,dk<R \meV{k) 




{2my 



F(sx,...,s k )Y[— ds 



( c fc) ( c i) 



i=l 



F(si,...,s k ) 



e' n 

b m , m£V(k) meV(k) 



E- 



Ll+r m I 



No ne o f the b m with m a singleton set are constrained by the divisibility conditions 
in ((^4), and other 6 m 's may also not be constrained when j s — j t and s 7^ t. Let 
Qj(k) denote the collect ion of m G 'P(fc) for which b m is not constrained by the 
divisibility conditions in ((^4); we will specify Qj(k) precisely below. We now write 
F(s±, . . . , s/c) as an Euler product 



(6.8) 



F(si,...,s k ) = H(l+ ^IT^r + /jfe s i' s 2,...,5 fe )) 



meQj(k) 



The local factor fj(p; s%, S2, ■ ■ ■ , s k ) — if p > max \j s — j t \, so we factor out the 
zeta factors corresponding to the unconstrained variables and write 



(6.9) F( S1 , . . . ,s k ) = [] C(l + r m ) { -^* m J hj( Sl , 

\m6Sj(fc) 

where we see hj(s±, S2, ■ ■ ■ s k ) is analytic in the region 



(6.10) 



C 



Writ,- 



1 < j < k. 



Our analysis is complicated by the different cases of non-distinctness in the 
numbers ji,j2,---,jk- Suppose there are r distinct values among the ji's. By 
reordering we may take these r distinct values to be ji, ]2, ■ ■ ■ , j r which occur with 
multiplicity 01, d2i • • • j o, r . We then have among the k original j's the relations, for 
1< i < r, 



(6.11) ji = j S2{l) = j S3(l) 



■ j Sa .(i), t < sz(i) < s 3 (i) < ■ ■ ■ < s ai(l) < k. 
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We also for notational purposes let 
(6.12) si(t)=i, 



l<i<r. 



Our divisibility conditions T>{j) now reduce to the divisibility of the numbers j u ~jv 
when 1 < u < v < r, which will be divisible by all the 6m's with m any element of 
the set 

(6-13) -P' M (k)={J (J V {Sm(u) , Sn(v)} (k). 

n—l m—1 

Hence T>(j) is given by the equations 

1 < u < v < r. 



(6.14) 



n ^ 

m£l" { „ : ,)(l) 



The indexing set Qj (k) for the unconstrained variables is now given by 
(6.15) 

r 

Qj( k ) = IJ 2*' where Qi = i m ^ ^ : m C { s l(*)> s 2(«): s 3(«)> ■ ■ ■ , S ai (i)}} , 
i=l 

where we see the Qi are disjoint from each other. As we saw in the cases of k — 2 



and k = 3 we will approximate F(s\, . . . , S&) in (6.9) by replacing each zeta factor 
with its first Taylor or Laurent term except when m is a singleton that does not 
occur as a subset of any other ml G Qj(k), in which case we retain the zeta factor. 
This latter situation occurs only when at = 1 and in this case Qi = {{?}}■ By 
reordering we can assume the multiplicities di form a non-decreasing sequence, and 
we define u> to be the number of dj's for which Oj = 1, and take lo — if there are 
no such <n's. Thus 

(6.16) a% = a,2 = • ■ ■ = a u = 1, 2 < a w+ i < a w+2 < . . . < a r . 

With this preparation, we now define fc(si, S2, ■ ■ ■ , s&) by 



(6.17) F( Sl) ...,s k ) = k( Sl , 



>(n«^ 



meC 3 -(fc) 



We thus see k(si, . . . , Sk) is analytic in the region ( 6.1C| ) except possibly when 
T m = 0, and therefore by ( |6.5|) 

(6.18) 



(2ttz) 



/••• / / Hs lt ..., 8k )( n ( 



(cfc) (c 2 ) (ci) 



A fl s ' d Si 



n ^j ds * 

i—uj+l 



As before we see that in the region (6.10) if |r m | 3> n»£m togHin+a) ^ or m e 
then 

(6.19) F(s u s 2 ,...,s k ) « fc exp( J D(log2||j||) 5 ) J] log B (| S4 | + 2) 

Ki<k 
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where the j dependence may be estimated as in (5.6). We first move the contours 
(ci), (02), (c w ) to C where we encounter only simple poles at s, = 0. For 
the remaining contours we proceed to move each (cj) to Cj and encounter various 
poles which generate new terms which need to be integrated over the remaining 
contours. In these new terms we continue to move each (cj) either to the left to Cj 
or to the right to the reflected contour Cj* according to whether the factor R aSj has 
a > or a < 0. This process eventually leads to a series of residues and integrals 
which are over contours Cj or Cj* . These integrals are all estimated as before. We 



arrange the contours so that they are well spaced where the conditions for (5.19) 
are satisfied, and apply (2.27) to see these integrals are, for log ||2j|| <C logi?, 



We conclude that 



Vj (fc) 



Tk(j) 



V, {j) log 4 R + O k {e 



i=0 



Here T>i(j) <Cfc exp(D(log 2| \j\\) s ) since these terms are obtained from derivatives of 
hj which will satisfy the same type of bound as in ( |5.6| ) . The leading residue terms 
will be those terms that contain the factor fc(0, 0, ... ,0), and hence are obtained 



by hj(0, . . . , 0) in the integral in ( |0.18 ) 



by replacing the factor fc(si, S2, ■ ■ ■ , Sk] 

Because the Qi's are disjoint, the new integrand without fe(si, S2, ■ ■ ■ , Sfc) may be 
broken into a product of terms with disjoint variables which may be integrated 
separately, and thus the main term is 



%(o,...,o) n 



1 



l +1 l (2^)"' 



n^) M)#w,+1 )n^ 

meQ, ■ ' 



(6.20) 



(Ca,) ( C1 ) 

1 f R s 



2iri J s 2 C(l +s) 



ds 



On moving (c) to C we see by (2.27) that 
1 f R s 



2ni J s 2 C(l + s) 

(c) 



ds = l + 0(e 



-c'\/Iog R 



and therefore these factors may be ignored in (6.2C) with an error of 0(e~ c Vi°g-R)_ 
As before we make a change of variable to remove the R dependence in the first 
integral but retain the Sj's as our variables in the r m 's through the replacement 
Sj — > Sj/ log R and (cj) — > ( . j R ). We thus have the main term above is 



(6.21) 
where 

(6.22) 



C fc (j> j (0,...,0)log^W R, 

r 

Cfc (■?)= n 



i=u)-\-l 
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(6.23) C k =C k (k) = 



1 



[2mY 



n ^) ( - i)#m+1 )n^^- 



me'P(fc) 



(<*) (ci) 

The power Vj{k) of log R obtained by the change of variables is equal to the number 
of r m in the denominator of the integrand minus the number of r m in the numerator 
plus k — to from the product of R Si / si 1 . Thus 



(6.24) 



u j (k) = k-w+ J2 (~ r - 

me Qj (fc) 
= fc+ ^ (-1)#»\ 



#m 



To evaluate the sum above, we use ( p. 15 ) and count the elements of Qj(k) 
according to their cardinality. Thus 



E 



(6.25) 



£(.-(1 

r 



+ ••• + (-1) 



ft, 



1=1 
-r, 



and therefore we have 

(6.26) Vj {k) = k-r. 

Next, from ( |6.8| ) and ( |6.9| ) we have, using ( |6.25 ) 

E 



ft i (o,o,...,o)=n(i 



1 ^meOj(fc) 



fi + - E (-i) #m + /i(p;o,o,. 



Letting 
(6.27) 

we will now prove that 
(6.28) 



l-- + / j (p;0,0,...,0) 
P 



f = (jl, J2, ■ ■• jr), 

r - v p {j') 



/i(p;0,0,...,0) 



and thus by (1.5) we conclude 

(6.29) ft i (0,0,...,0) = 6(j / ). 

From @, and deTTj) we see that 

/i(p;si,s 2 ,-..,sfc) = E E 

l<u<D<r m£ , P'j„ i „ ) (i) 

p|j'«— iti 



(-!)#'" 
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and hence 



/;(p;0,...,0) = 1 £ E 



Ku<«<rtBep' {Uit j(fc) 

pIj'd— j'« 



Thus, the proof of ( |6.28 ) reduces to proving 



E E (-i) #m =r--^on 



If v p {j') = q, then the r distinct numbers j\, j'2, ■ • ■ , jV must fall into q residue classes 
congruent to numbers r% , r-x , . . . , r q modulo p, say with multiplicities uii , m,i , . . . , m q , 
with rrii > 1 and mi + to 2 + ■ ■ ■ + Ti ? = r. Let 

Mi = {n : ] n = r. t (mod p), 1 < n < r}, 



so that #Mi = m,. Next, recalling ( 3.15 ), let 

Mi = {si(n),s 2 (n), . . . ,s „(n) MJ 

= U U s »( n )' 

and 



={m^|:mC Mi, #m > 2}, 



so that 



#M, ; = ^ a„. 

The condition p\j v — j u will only hold if j v and j u are in the same residue class, 
and hence if and only if 

{u,v}e\Jv(i). 



i=l 



Hence we conclude 



l<?z<t;<r 



since "P(i) also contains all the unconstrained sets of variables with > 2 elements 
which are not included in 7-"{ u ,v}{k) and Qj also includes the singleton sets which 
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are not in V{i). We conclude by (3.25) 



E E (-!) #m =E E (-i) #m - E (-i^-E 1 



l<u<v<r meV'{ u , v y(k) 

P\ji>-ju 



i=1 meP(i) 

EE (*f 

i=l j=2 v J 

E ( - 1 + #m 



(-iy +r- k 



- k 



= -q + E E a n + r - k 

= r — q 

= r- v p (j'), 



which proves ( |6.28 ). 

Combining our results we have shown that 



fc-t 



(6.30) T k (j) = C k (j)6(j') log k - r R + J2 Viti) ^g k ~ M R + O k ( e - c '^^) 



Returnin g to t he notation of (1.2), we see in ( |6 . 2| ) that S(j) = Sk(N,j', a) wit h j' 
given by ( |6.27|) . Thus we have proved the following refined version of Theorem |l.l[ 

Theorem 6.1. Given k > 1, let j = (ji , j2 , ■ ■ ■ ,jr) and a — (ai,a,2, . . .a r ), where 
the ji's are distinct integers, and a,; > 1 with J2i=i a i — k- Assume log2||j|| <S 
logi?. Then we have 



(6.31) 



k—r 



S k (N,j,a) = C k (a)&(j)N(logR) k - r + ^ V l (j)Nlog k - r - 1 R 



+ O k (\\f\\Ne- c '^^) + 0(R k ), 



where V t (j) < exp( J D(log2||j||) 5 ). 



The constants C k (a) satisfy the relation (1.7) because we may extend the product 
in (6.22) to include the terms with 1 < i < li since for these terms C ai = C\ = 1. 



7. The constants C k 

We now discuss the constants C k that arise in the fc-correlation results. We will 
show the integrals defining these constants converge absolutely and vanish as the 
contours are moved to the left to — oo. New terms generated from the residues can 
be handled similarly by moving the contours appropriately to the left or right to 
infinity. Therefore the C k 's can be calculated by a residue calculation. We have 
already computed C2 and C3 before, but we will recompute them directly and also 
compute C4. 
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Consider first 

(7.1) C 2 =C 2 {2) = —— 2 / / d Sl ds 2 , 

{2-Kiy J J SlS 2 {Sl+S 2 ) 

(C2) (Cl) 

for c 2 > ci > 0. Letting 

(7.2) J(cri,<7 2 ) = ^2^)2 y y m 2(si, S2) dsi rfs 2 , 

("•2) (cti) 

where 

(7.3) m 2 (si,s 2 ) = < ; r, 

SlS 2 (Sl + s 2 ) 

we may rewrite ( |7.l| ) as 

C 2 = I(<7i,a 2 ), a 2 > <7i > 0. 

If 

(7.4) N>1, N>1, |si + s 2 |>l, 
then we have the bound 

„cri+cr 2 

(7.5) m 2 (si,s 2 )< 



(l + |*i|)(l + N)(l + |ti + *a|)' 
Suppose the contours (oi) and (<7 2 ) are arranged so that the conditions in ( [7.4[ ) 
hold; this can always be done in what follows. Then by symmetry we may estimate 
I(&if&2) over the restricted integration range |si| < \s 2 \ and t 2 > 0, and hence 



i+* 2 vy-* 2 (i + N)(i+*i+*2) 

Jo (i + ^J 2 

< e CTl+<T2 . 

From this estimate we see C 2 is well-defined, and further that 
(7.7) lim l(a 1 ,a 2 )=0, j = 1,2. 

(Tj — * — OO 

Since m 2 (si, s 2 ) — » if \ti\ — > 00 or |i 2 | — > 00, we can move the contour integrals 
defining C 2 to the left and C 2 can be evaluated from the residues encountered. 

Moving first c\ — > —00, we encounter residues at the simple poles s\ = and 
si = — s 2 , and hence 

1 f f e &2 1 \ 

0=) 

Next moving c 2 — > —00, we encounter a pole at s 2 — with residue equal to 1. 
Hence 

(7.9) C 2 (2) = 1. 
Turning next to C 3 , we have for 03 > <r 2 > <7i > 

(7.10) C3=C 3 (3) = 7(a 1 ,a 2) (T3) 
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I {ax, 0-2,(73) := / 2 ^\3 y y J m3{si,S2,s 3 )ds 1 ds 2 ds 3 

TO 3 (S1,S2,S 3 ) 



(0-3) (0-2) (<ri) 

(si + s 2 + s 3 )e Sl+S2 + S3 



Sis 2 s 3 (si + s 2 ){s 1 + s 3 )(s 2 + s 3 ) ' 



N > 1, 1 < i < 3, I + sj\ > 1, 1 < i < j < 3, 



where 
(7.11) 

(7.12) 
If 

(7.13) 
then 

and hence 
(7.14) 

|si||s 3 ||si + s 2 \\s 2 + s 3 | 

By symmetry it sufheies to estimate / over the region \si\ < \s 2 \ < |s 3 |, and hence 
as in (7.6), assuming (7.12) holds, 

(7.15) 



si + s 2 + s 3 



1 



s 2 (si + s 3 ) s 2 Si + S 3 

m 3 (si,s 2 ,s 3 ) « e - 1 +^+^ 



«1, 



1 



/(<ti,<7 2 ,<7 3 ) < e 



(^3) 



< e 



*3_ 

M 



1 



S2 + S 3 



1 



Si Si + s 2 



O2) 

|s 2 |<s 3 



log(2+M) 



(o-i) 

KI<M 

dU 



\S 2 \\S 2 + S3 



(*s) 



(<*0 

l*2|<s 3 



< e 



f 1 + C"2 + C"3 



log 2 (2+|t 3 |) 



S3 



(^3) 



< e" 

Thus C 3 is well defined and 

(7.16) lim l{a 1 ,a 2 ,a 3 ) = 0, j = 1,2,3. 



(Tn — ► — OO 



Since m(si, s 2 , s 3 ) vanishes as \ si\ — ► 00 for any 1 < i < 3, we can move our contours 
to the left and evaluate C 3 from the residues encountered. 

We now evaluate C 3 by making a simplification suggested by D. Farmer which 
simplifies the calculation of the residues in this case. Observe that if we break the 
integral defining C 3 into three pieces corresponding to each term in the numerator, 
then by symmetry these terms are equal to each other, and therefore 



(7.17) C 3 



{2mf 



a Si+S 2 +S3 



■ dsi ds 2 ds 3 



s 2 s 3 (si + s 2 )(si + s 3 )(s 2 + S3) 

(c 3 ) (c 2 ) (ci) 

We choose Cj — j for 1 < j < 3. Moving first c\ — * —00, we encounter the simple 
poles si = — s 2 , and si = —S3 and obtain 



C 3 



(2-Ki) 2 



(c 3 ) (c 2 ) 



s 2 s 3 (s 2 + s 3 )(s 2 - S3) 



• ds 2 ds 3 . 
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As a function of s 2 the integrand has simple poles at s 2 = 0, s 2 = S3, and s 2 = —S3. 
On moving the contour c 2 — * —00 we do not encounter the pole at s 2 — S3, since 
S3 is on (03) = (3) and is to the right of (C2) = (2). The residue at S2 — is 

e S3 - 1 



S3 3 



and the residue at S2 — — S3 is 



therefore we have 

(oa) (c 3 ) 

In the first integral above we move C3 — > — 00 and encounter a triple pole at S3 = 
with residue equal to j. The second integral is equal to zero since we encounter no 
singularities when we move 03—5-00 where the integral vanishes. We conclude 

(7.18) C 3 = \. 

We now turn to the general case, and have for u\ , o 2 , ■ • ■ , (Jk > 

(7.19) C fe =4(0-1,02, 



where 

(7.20) 4(0-1,(72, ... ,0-fc) = t^twT I ■■ I m k (s 1 ,S2,-.-,Sk)ds 1 ds2---ds k , 



1 

(2tt7)' c 



(7.21) m,( S i, S2 ,... Sfe )= ( [] (r m )(- 1 ) #m+1 ) JJ ^ 

/ Si 



meV(k) i=l 



and r m is defined in (6.7). We will prove that 4 is well defined and converges by 
proving that provided |r m | 3> 1 for m € 'P(fc) then 

(7.22) 4(oi, 02, e "i+«a+-+^. 

We have already proved this for fc = 2 and fc = 3 and we will proceed induc- 
tively Define a function gk{s\, S2, . . . , Sfc) by first taking 52(si, S2) = 1, and define 
53(si,s 2 ,s 3 ) by 

gSl+S2+S3 

m 3 (s 1 ,s 2 ,s 3 ) = — ■ ■ r— 3 3 (si,S2,s 3 ) 

si(si + s 2 )(s 2 + s 3 )s 3 

so that 

, v si + s 2 + s 3 

53(si,s 2 ,s 3 ) = — ■ r-. 

s 2 (si + s 3 ) 

We next define gi(si, s 2 , S3, S4) by 

gSl+S2+S3+S 4 

m 4 (si, s 2 , s 3 , s 4 ) = — ■ ■ ■ r — ga(si, s 2 , s 3 )g 4 (si, s 2 , s 3 , s 4 ) 

si(si + s 2 )(s 2 + s 3 )(s 3 + s 4 )s 4 

so that 

(si + s 2 + s 4 )(si + S3 + s 4 )(s 2 + S3 + s 4 ) 



54(Sl,S 2 ,S3,S 4 ) 



s 3 (si + s 4 )(s 2 + s 4 )(si + s 2 + s 3 + s 4 ) 
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Wc define g k in general by 



fh-x 



(7.23) m fc (si,s 2 , ■ ■ ■ ,Sfe) = (n^* J 7777 ( II J. + s - +1 I Y[9i(si,a^,...,Si). 



i=3 



We will show below that, provided \r m \ ^> 1 for to £ V{k), 



(7.24) 



\\\9i{si, s 2 , ■ ■ ■ , s k )\ < fc 1. 



4=3 



Assuming this estimate for the moment, we see 
(7.25) 



As in ( |7.6| ) and (|7.15| ), provided |sj| 3> 1 and |s< + Si+i| >lwe see that 
(7.26) 



l«K i '(-> + ti).^. < log B+1 (2+|t i+1 |) 



M<|s i+1 | 



By symmetry we may estimate the integral in (7.25) over the region |si| < | S2 1 < 
• • ■ < \sk\, and hence we see on using the estimate ( 7.26]) re peatedly that the integral 
converges, which proves ( gjg ). It remains to prove ( [7.24| ). 

We have already seen that 53 <C 1 in the equation above (7.14). For g 4 we have 

Isi + s 3 + s 4 | |si + s 2 + S4IIS2 + s 3 + S4I 



(54(51,52,53,84)1 



|S 3 ||Sl + s 4 | 

1 1 



S3 Si + S 4 

« 1, 

since r m 3> 1. For the general case we claim 
(7.27) 

, , (si + s fc _i + sj,) _ 

5fe(Sl, 52, ■ ■ • , Sfe) = ; ; 



1 52 + S4||Sl + S 2 + S3 + S 4 | 
1 J £J_ 



S2+S4 



S2+S3 + S4 



Sfc-l(si + Sfe) 

(si + s 2 H h Sfe_ 2 + s fc )(s 2 + s 3 + 



Sfe, 



, (s 2 + s 3 H h s fe _ 2 + Sfc)(s a + s 2 + 

and gfc is a product of terms of the form 

(-!)*» 



Sfe, 



(-1)" 



(7.28) 



si + r/ m 

?7m 



si + rj m + Sj 



Vn 



E 



where to S V{k), to has between two and k — 3 terms, and fc £ to but 1 ^ to since 
we have separated out these terms above. For each to we may choose j arbitrarily 
except for the requirements that 2 < j < k and this factor has not already appeared 
in the product. For example, 

(si + s 2 + s 5 )(s 2 + s 4 + s s ) (si + s 3 + s 5 )(s 3 + s 4 + s s ) 



95 



(s 2 + s 5 )(si + s 2 + s 4 + s 5 ) (s 3 + s 5 )(si + s 3 + s 4 + s 5 ) 
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For k = 6 we can write ge as a product of 6 terms corresponding to the pairs (f] m ,j) 
given by (s 2 + S6,3), (s 3 + s 6 ,4), ((s 4 + s 6 ,2), (s 3 + s 5 + s 6 , 2), (s 2 + S5 + s 6 , 4), and 
(S4 + S5 + Se, 3). In the general case, we need to check that we will run through all 
the m £V(k) with this decomposition. Then, since 



(7.29) 



si + r)„ 



Vn 



Si + Tim + Sj 



1 + -23- 

»7m 



<fc 1 



(-1)* 



if |t to ] 3> 1, and similarly for the first and last individual terms in ( 7.27 ), we see 
that ( 7.24 ) follows. To establish that gk can be written as a product of terms of the 
form (7.28), we use a counting argument to verify that in constructing the product 
we use all of the remaining terms r m with m 6 V(k) in TO/-(si, s 2 , . . . ,s/c). We 
examine this by grouping according to #to. All of singleton elements have already 
appeared in ( |7.23| ). Note every m for which rj m occurs in g^ must contain k, because 
all the terms r TO without Sk have already appeared in some gi, 2 < i < k — 1, or the 
product in ([7.23|) . Further, 1 can not be in m because we have separated out terms 
with si in (7.28). Therefore, for rj m with #m = 2, we can take r\ m = Si + Sk with 
2 < i < k — 2 since Sk-i + Sfc already appears in the product in (7.22), and hence 
there are k — 3 such terms. All of #m = 2 terms are now accounted for in (7.23). 
For r\ m with #771 = 3, we can take i^^ 2 ) choices for rj m in the first denominator 
of ( [7.28 ) except we can not use any of the k — 3 factors r\ m + Sj we just obtained 
when #m = 2, and hence there are 



2 H fc - 3 )= 2 



such terms in (7.28), and this accounts for all the r m with #to = 3. In general 
we see that there will be (*Z^) terms 77 m with ^=771 = j occurring in (|7.28| ) for 
1 < j 5: ^ ~ 3. Since 7^777 = fc — 3 is the largest size set for which rj m occurs in 
( 7.28| ), we only need to check that the terms with j^m = k — 2, j^m = k — 1, and 
#777 = k are accounted for in the earlier terms and the last term in (7.27). There is 
only one element with #777, = fc, and it only occurs in the last term in (7.27). For 
the k — 1 elements with #777 = k — 1, k — 3 of them occur in the product formed 
from elements with 7^777 = k — 3 in ( |7.28[ ), and the remaining two terms are in the 
last term of (7.27). Finally, for the ( 2 ) elements with #777 = k — 2, we have (jjlj?) 
of them in the terms in ( 7.28 ) coming from terms rj m with #777 = k — 4 and 2(fc — 3) 
of them from rj m with #777 = /« — 3. Since 



k - 1 
2 



/c-3 
fc- 5 



2(fc-3) 



we see this leaves exactly the one term needed for the last term in ( 7.27 ) 
We conclude this section by evaluating C4. We have 



(7.30) 



C 4 = ^2^)4 J J J J m4 ( Sl ' S2 ' S 3,S4) ds 1 ds 2 ds 3 ds4, 



(04) (c 3 ) (c 2 ) (ci) 
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where C4 > C3 > C2 > ci > and 

TO 4 (si, S 2 , S3, Si) = 

(fi + s_2 + ggKgi + 52 + s 4 )(si + S3 + s 4 )(s 2 + s 3 + s 4 )e Sl+S2+;i3 + s - 
s 1 s 2 s 3 S4{si + s 2 )(si + s 3 )(si + s 4 )(s 2 + s 3 )(s 2 + s 4 )(s 3 + s 4 )(si + s 2 + s 3 + s 4 ) ' 

Moving ci to — 00 we encounter poles at si = 0, Si = — s 2 , Si = —S3, si = — s 4 , 
and si = — s 2 — S3 — s 4 . The residues at the poles si = — s 2 , Si = —S3, si = — s 4 
are symmetric with each other, but because the contours are not symmetric they 
do not contribute equally. 

First, the residue at s± = is equal to 



s 2 2 s 3 2 s 4 2 



and thus by (2.8) or directly we see immediately that this term is equal to 1. 
Next, the pole at si = — s 2 — S3 — s 4 gives a residue equal to 

1 



(s 2 + S3) 2 (s 2 + s 4 ) 2 (s 3 + s 4 ) 2 



Moving c 2 to —00 we pass poles at s 2 = —S3 and s 2 = — s 4 and the two residues 
here cancel each other out and we get a contribution of 0. 
For the pole at si = — s 2 we have a residue 

(s 2 - s 3 - 34) (s 2 + s 3 + s 4 ) e 33+3i 

s 2 2 (s 2 - s 3 ) (s 2 + s 3 ) (s 2 - s 4 ) (s 2 + s 4 ) (s 3 + s 4 ) 2 

Moving next c 2 to —00 we pass poles at s 2 = 0, s 2 = —S3 and s 2 = — s 4 . The 
residue here at s 2 = contributes 0. The residue at s 2 = —S3 is 

s 4 (2s 3 + s 4 )e S3+s * 
2s 3 3 (s 3 - s 4 )(s 3 + s 4 ) 3 ' 

which on moving C3 to the left passes poles at S3 = and S3 = — s 4 ; on moving c 4 
to the left the first pole gives a contribution of — i and the second gives 0. In a 
similar way we find the pole at s 2 = — s 4 contributes 0. 

Returning to the poles at s 2 = — S3 and s 2 = — s 4 , we will see presently they 
contribute 0. Both poles can be handled similarly, so we consider only s 2 = —S3. 
We move c 2 to the left and encounter poles at s 2 = —S3 and s 2 = — s 4 . The pole 
at s 2 = — S3 gives a residue 

s 4 (-2s 3 + 3 4 ) e - a3+fl * 
2s 3 3 (-s 3 + s 4 ) 3 (s 3 + Si) ' 
We move C3 to the right to make the exponential factor vanish, and encounter only 
the pole at S3 = s 4 which contributes the negative of the residue 

3 + 6s 4 + 2s 4 2 
16s? 

which on moving c 4 to the left makes a contribution of 0. In the same fashion the 
pole s 2 = — s 4 makes no contribution. 

The poles at Si = —S3 and si = — s 4 can be handled as above, and one finds 
that they contribute 0. Hence 

(7.31) C 4 = \. 
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8. A Generalized Result 

To obtain the best numerical results in applications one needs the truncation level 
R as large as possible, and this can be optimized by using different truncations Ri 
in the correlations. Thus, in place of (6.1) we consider the sum 

N 

(8.1) S{j,R) = ^A fll (n + j 1 )A i i 2 (n + j 2 )...A flfc (n + j fe ), 

n=l 

and see as before 

S ^ =N 22 [ dl , d2 ,...,d k ] +0( R R ^ 

di<Ri, l<i<k ' 1 i=l 

(8.2) (d r ,ds)\js~j r ,l<r<s<k 

k 

= NT k (j ) R) + 0([[R i ). 

i=i 

We now let 

(8.3) Ri = N 0i , l<i<k, 
and let — (6i, 02, ... ,6k)- Then ( |6.5|) becomes 

(8.4) r fc(j , jR ) = __ I y ... J F ( Sl ,..., Sfc) p__ ds . 

( c fc) ( c i) I_1 

The analysis is now identical to before with the additional constants 9i only effecting 
the constants obtained in the main terms. Suppose as before j = (ji,j 2 , • • ■ jk) 
has r distinct values which we may take to be the first r components of j with 
multiplicities a = (ai, a 2 , ■ • ■ , a r ). Then the truncations = (9\, 6 2 , ■ ■ ■ 9 k ) will be 
partitioned into r sets. We denote by 9i(v), 1 < v < ai the 0's associated with ji, 
and write 

= [0 X ,0 2 , ■■■ , r ), t = (0,(1), 9,(2), ■■■ , Ofa)). 

Let 



N r / m 

.5) S k (N,j,a,0) = J2T[ I1 A ^ 

n— 1 i— 1 1 



(»)(n + ji) j , 

Theorem 8.1. With the above notation, we have for < Oi < 1, and log2||ji|| <C 
log AT 

(8.6) S k (N,j,a,0) = (C k (a,0)e(j) + o k , e (l))N(logN) k - r + 0(N e i +e >+--- +9k ), 
where 

r 

(8.7) C k (a,O)=\\C ai {0i), 

i=l 

and 

(8.8) c fcW = * /•-/( n ^) ( -^ +i )n^- 
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We see that in the special case 9 = (6,6,- ■ -6) and Ri = N e that 

(8.9) c k (o) = e k - l c k . 

In the general case we can compute these constants by the same residue calculations 
used in the last section. The results depend on the relative sizes for the 0j's, and 
we find using the ordering Oi > 9 2 > ■ ■ ■ > 9 k that C\ = 1, and 

(8.10) c 2 {e) = e 2 = mm{e 1 ,e 2 ), 

(8.11) C 3 (6) = 2 03 - \(e 2 + 3 - 0i) 2 [0 2 + 3 > 0i], 
and, letting A 3 = 9 2 + 3 - 6 1 and A 4 = 6\ - 6 2 - 6 3 + 6 4 , 

(8.12) c 4 (6>) = e 2 e 3 e A - U a a 3 2 \a 3 > o] - ^a a (a 4 2 + 6A 3 a 4 + 4A 3 2 ) \a a > o]. 
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